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INTRODUCTION 


^ The  Least  Mean  Square  (LMS)  algorithm  is  a method  of  extracting  a 

l9" 

target  signal  in  the  presence  of  interfering  noise  sources.  WidroW 
£- 

and  Griffiths  developed  the  theory  and  applied  simulated  tests  using 
random  numbers.  These  tests  indicate  that  the  method  is  feasible  for 


Laboratory  conditions.  The  nurrose  of  the  sea  test  is  to  evaluate 


the  method  for  actual  ocean  conditions.  In  this  section,  the  exrerL 


raent  of  adaptive  beamforming  is  described.  This  includes  a brief  des 
crintion  of  adaptive  beamformingjthe  computer  "rogram^and  the  results 
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EXPERIMENT  ' 

The  geometry  of  the  experiment  is  given  in  Figure  1.  A four  element 
vertical  array  is  suspended  from  the  FLIP.  The  depths  of  the  elements 
are  2,500,  2,521,  2,563,  and  2,601  ft.  The  USNS  SANDS  (AGOR-6)  is  at  a 
distance  of  10,000  ft  from  the  FLIP.  A 400  Hz  projector  is  suspended 
from  the  USNS  SANDS,  at  a depth  of  1000  ft. 

The  Honeywell  projector  is  driven  by  a pseudo-random  noise  (PRN) 
signal  centered  at  500  Hz  in  a 200  Hz  band.  The  HX-90  has  a maximum 
source  level  of  100.6  db  (/TC  lL>bar  at  1 yd).  Driving  the  projector 
with  a PRN  signal  reduces  the  level  to  91.1  db  and  at  a range  of  10,000 
ft  the  expected  acoustic  levol  at  the  receiving  hydrophones* is  approx- 
imately +19  db.  The  frequency  response  of  the  HX-90  is  shown  in 
Figurd  2. 

A block  diagram  of  the  transmitting  system  is  shown  in  Figure  3. 

A PRN  signal  that  can  either  be  continuous  or  gated  is  filtered  through 
a Butte rworth  filter  that  has  a 200  Hz  bandwidth  centered  at  500  Hz. 

The  signal  is  then  amplified  and  transmitted.  The  CML  amplifier  has  a 
maximum  output  power  of  5 KV  and  with  a 40$  projector  efficiency,  the 
electrical  power  is  more  than  adequate  to  drive  the  projector  at  full 
power  if  required. 

A block  diagram  of  the  receiving  equipment  is  shown  in  Figure  4> 

The  signal  received  at  the  hydrophone  is  amplified  and  the  information 
is  telemetered  back  to  the  ship  where  it  is  filtered  through  identical 
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Buttervorth  filters  and  recorded  on  magnetic  tape  in  analog  form.  The 
analog  signals  also  go  to  a 12  bit  A/D  converter  and  these  are  processed 
on  the  1230  Univac  computer.  The  digitized  signals  are  also  recorded  on 
Uni vac  1540  magnetic  tape. 

In  order  that  the  projector  be  in  the  Fraunhofer  region  of  the  array, 
the  distance  D between  projector  and  array  must  be  such  that, 


0 £ ^ 


where  is  the  length  of  the  array  and  A is  the  wavelength. 

For  500  cps,A~  10  ft.  Letting  Jl  = 100  ft.,  then  D > 5000  ft. 

Since  D = 10,000  ft.,  the  signal  from  the  projector  may  be  approximated 
by  a plane  wave  at  the  array.  The  beamwidth  (bw)  of  a linear  array  for 
a conventional  pattern  is 


>SS: 

Thus,  for  our  experiment,  the  bw  ^ 6 degrees. 


Adaptive  Beamforming 

1 2 

We  will  follow  the  theory  developed  by  Widroa  and  Griffiths  on  the 
subject  of  adaptive  beaoforming.  A tapped  delay  line  is  attached  to 
each  element  of  a K element  array I There  are  L weights  in  each  tapped 
delay  line.  (See  Figure  5)«  The  purpose  of  the  algorithm  is  to  adjust 
the  .weights  such  that  there  is  a least  mean  square  (LMS)  difference 
between  the  output  y and  the  expected  signal.  The  problem  was  first 
solved  by  Wiener  and  the  final  answer  was  in  terms  of  an  inverse  matrix. 
Widrow  and  Griffiths  obtained  an  approximate  iterative  solution.  The 
main  attractive  features  of  the  solution  is  that  the  noise  field  need 


I 

j 
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In  our  simulation,  the  power  spectrum  of  the  signal  is  assumed  to  be 
known.  In  obtaining  the  LMS  difference,  the  noise  tends  to  be  minimized 
and  the  pattern  of  the  array  tends  to  have  notches  where  the  noise  inter- 
ferences are  strongest. 

Based  on  the  thermal  structure  of  the  area  and  a ray  tracing  program. 


the  angle  from  the  array  to  the  propagation  path  of  the  plane  wave  from 
the  projector  is  found  to  be  80°.  This  is  the  direct  path.  In  the 


computer  program,  the  array  is  steered  broadside,  so  that  a signal 
coming  from  this  direction  is  enhanced.  A signal  is  not  simulated  in 
this  experiment.  A plane  wave  striking  the  array  in  any  direction 
other  than  broadside  is  called  an  interfering  noise  in  this  article. 
Thus  the  plane  wave  from  the  projector  is  an  interfering  noise  and 
one  of  the  purposes  of  this  test  is  the  nulling  out  of  this  noise.  It 
should  be  noted  that  the  interfering  directions  are  not  inputs  to  the 
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second  taps  of  the  k receivers  etc.  Thus  at  any  iteration  there  is  a 
set  of  Z\.  » where  i = 1,  2,  ...  KL.  Let  the  center  point  of  each 

tapped  delay  line  be  the  reference  point  for  the  correlation  values.  * 

. 

For  an  even  number  of  taps,  the  correlation  of  the  signal  between  the 
adjacent  tap  and  the  reference  point  is  (4)  * where  4 is  the  delay 
shown  in  Figure  5»  Using  the  symmetry  properties  of  time,  the  auto- 
correlation  of  the  input  to  the  first  tap  with  that  of  the  reference 
point  is  ^ ^ Thus  a table  of  KL  values  of  computed  and 

stored  in  the  computer. 

The  computer  program  used  in  this  test  includes  two  equations.  The 

i 

first  equation  updates  the  weights  for  each  iteration.  After  each 
1000  iterations,  the  beam  pattern  subroutine  is  called  and  these  weights 
are  used  to  compute  the  beam  pattern.  The  iteration  equation  is 

CO  fUO-y  Sjntej], 

ce  U (1) 

y - wL^xii) 

L-! 

where 

x(i)  is  the  input  at  the  ith  tap 
- w(i)  is  the  weight  at  the  i^  tap 
y is  the  output  of  the  array 

Ps(i)  is  the  signed  correlation  between  the  i^  tap  and  the 

reference  point 

, X<0 
' - + I ) X > o 

\j  is  L,  r~  5 
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Equation  (l)  by  itself  is  not  a difficult  equation  to  program.  The 
inclusion  of  fixed  point  arithmetic  and  subroutines  for  real  time  pro- 
cessing increase  the  difficulties. 

For  an  array  steered  in  the  0^  direction  and  for  a plane  ware 
striking  at  the  & direction,  the  output  of  the  array  is 

■>*.■=/  L~° 

where  'j-  is  the  delay  time  to  reach  the element. 

f/*s 

The  power  output  of  the  array  is 

U,(e)=  - 27 

>yzi  i-i  £-• 

+■'*&) -7X.(e.,)  +TXce)-^w7  ti) 

where  Rtf') is  the  autocorrelation  function  of  the  stationary  random 
process  idt).  There  are  terms  to  compute  for  each  value  of  Q . 

The  time  required  to  compute  equation  (2)  can  easily  be  in  terms  of 
minutes  for  the  case  of  180  values  of  Q . In  order  to  overcome  this 
difficulty,  terms  of  identical  arguments  are  regrouped.  A complete 
description  of  this  regrouping  and  taking  advantage  of  symmetry  proper- 
ties is  given  in  reference  (3).  The  resulting  equation  is 

ft  i C «-*-»,  ) 

1*~o  -»*•*/ 

^ 2.  Z O)  r Oj  s ) (3) 

<?-S|  '**1*1+1  ' 

4-0.  t"' 2 2 cC—,nA)'r<:~,,nJs  0V,  6) 

P>-  I /**)*/  n~  I 

W-W  Gp-‘j  4=0 

■'  *■/  f-  > ° 
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Whereas  equation  (2)  requires  terms,  equation  (3)  requires 
K4- 1")-  ('£21^')  terms  of  ~r~C'r'i  nt  ‘ ©<*,  &)■ 

Thus,  we  compute  )2.iy  terms  for  each  $,  instead  of  1600  terms.  This  is 
the  minimum  number  of  values  of  T required  for  an  arbitrary  array 
geometry.  Although  equation  (3)  looks  more  formidable  than  equation  (2), 
the  programming  difficulty  was  on  the  same  order.  The  subroutines  in 
CS-1  language  are  given  in  reference  (4) , for  equations  (1)  and  (3). 

The  inputs  at  each  receiver  is  sampled  at  1,024  samples/sec.  The 
total  sampling  time  is  4 secs.  Thus  the  total  number  of  samples  for  the 
4 receivers  is  1024  x 4 x 4 = 16,384  samples.  The  projector  emitted 
the  PRN  noise  for  a period  of  6 secs.  The  computer  obtained  4 secs 
of  this  digitised  data  and  this  data  was  stored  in  memory.  The  A-D 
conversion  is  12  bits.  Each  word  of  memory  can  store  30  bits.  Since 
half  words  may  be  used  in  the  Univac  1230,  the  total  input  data  is  stored 
in  8,192  words  of  memory.  The  1230  has  32,000  words  of  memory,  so  that 
there  is  ample  storage  for  the  program,  after  the  data  is  stored.  In 
the  program,  the  number  of  taps  in  each  delay  line  is  10,  i.e.  L = 10. 
Results 

The  computer  program  was  first  simulated  on  ambient  noise  alone,  i.e. 
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the  projector,  called  the  interfering  noise,  was  not  turned  on.  The 
ratio  of  interfering  noise  to  ambient  noise  was  lOdb.  The  results  of 
this  experinent  for  one  run,  is  given  in  Figures  (6)  and  (7).  In  this 
case  ^ was  set  to  .005.  The  pattern  was  drawn  at  the  end  of  each  1000 
iterations.  There  is  no  strong  interfering  plane  wave  and  no  large 
notch  is  obtained.  All  patterns  are  siailar  to  one  another.  For  the 
next  case,  the  projector  was  turned  on.  The  results  of  this  case  are 
shown  in  Figure  (8)  and  (9).  For  this  run,  \J  - .01.  At  the  end  of 
1000  iterations,  there  is  no  definite  notch.  At  the  end  of  2000 
iterations  there  is  a notch  of  14«5  db  at  82°  and  at  the  end  of  3000 
iterations  it  is  17^  db  at  82°.  At  the  end  of  4000  iterations  it  is  15^ 
db  at  83°.  Each  run  for  the  case  of  the  turned  on  projector  resulted  in 
notches  siailar  to  the  case  shown. 

A notch  corresponding  to  the  surface-reflected  path  was  not  obtained, 

probably  because  the  test  was  conducted  during  very  rough  weather  and 

the  surface  reflected  signal  may  have  been  severely  attenuated. 

* « 

Discussion  and  Conclusions 

The  primary  conclusion  is  that  adaptive  beamforming  does  work  in  an 
ocean  environment.  The  conditions  in  the  medium  during  4 secs  of 
adaptation  are  evidently  stable  enough  for  the  adaptive  iterative  pro- 
cess  to  converge.  A notch  of  about  17  db,  an  improvement  of  approxl- 
nately  7 db  over  that  obtained  in  the  absence  of  the  interfering  source 
is  good,  realising  that  there  were  only  4 receivers  in  this  experiment. 

Since  there  was  only  1,024  samples/sec  and  since  the  highest  fre- 
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quency  passed  through  the  filters  was  600  Hz,  there  was  aliasing  and  the 
distortion  of  the  power  spectrum  and  in  turn  the  correlation  values. 

The  sampling  rate  should  have  been  greater  than  or  equal  to  1200  samples. 
Of  course,  we  could  have  sampled  at  twice  the  bandwidth  to  obtain  all 
the  information  in  the  input,  but  this  would  require  changing  the 
algorithm.  Since  there  was  distortion  in  the  power  spectrum,  then  the 
notch  obtained  was  not  at  a minimum.  Future  work  is  planned  to  find 
the  numerical  effect  of  aliasing  in  this  case,  on  the  notch.  It  is 
also  planned  to  write  a minimum  variance  distortionless  look  (MVDL) 
on  the  data.  This  will  result  in  a lower  signal-to-noise  gain  but  no 
distortion  in  the  signal.  It  is  also  planned  to  carry  out  the  adaptive 


beamforming  method  on  other  fixed  arrays  under  different  ocean  conditions. 


BENJAMIN  F.  CRON 


9 


NUSC-NL  Tech  Memo 
No.  2211-178-70 


REFERENCES 

1.  B.  Widrow,  P.  Mantey,  L.  Griffiths,  and  B.  Goode,  "Adaptive  Antenna 
Systems,"  Proc.  IEEE,  55 » 12,  Dec  1967 

2.  L.  Griffiths,  "Signal  Extraction  Using  Real-Time  Adaptation  of  a 
Linear  Multichannel  Filter,  "Technical  Report  #6788-1,  Stanford 
University,  Feb  1968 

3.  A.  Nuttall,  D.  Hyde,  B.  Cron,  "Efficient  Computation  of  Beam  Patterns 
for  Arrays  with  Tapped  Delay-Line  Combiners,"  NUSL  Tech  Memo  2020-88-69, 
May  1969 

4.  B.  F.  Cron,  and  G.  Botseas,  "A  Scheduled  at-Sea  Simulation  of 
Adaptive  Beamforming,"  NUSL  Tech  Memo  No.  2211-162-69,  Sep  69 


NUSC-NL  Tech  Memo 
No.  2211-178-70 


PROJECTOR  NOISE 
V 500  Hi 
tW  • 200  Hi 


3000  ITERATIONS 


FIGURE  9 


NUSC-NL  Tech  Memo 
No.  2211-178-70 


